ANALYSING THE ELASTICITY DIFFERENCE TENSOR OF 
GENERAL RELATIVITY 
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Abstract. The elasticity difference tensor, used in [T] to describe elasticity 
properties of a continuous medium rilling a space-time, is here analysed from 
the point of view of the space-time connection. Principal directions associated 
with this tensor are compared with eigendirections of the material metric. 
Examples concerning spherically symmetric and axially symmetric space-times 
are then presented. 



' 1. Introduction 



In recent years there has been a growing interest in the theory of general relativis- 
tic elasticity. Based on the classical Newtonian elasticity theory going back to the 
17th century and Hooke's law, some authors began to adapt the theory of elasticity 
to the relativity due to the necessity to study many astrophysical problems as the 
interaction between the gravitational field and an elastic solid body in the descrip- 
tion of stellar matter, as well as to understand the interaction of gravitational waves 
and gravitational radiation and to study deformations of neutron star crusts. One 
■ of the first elastic phenomenon considered in the relativistic context was Weber's 

observation of the elastic response of an aluminium cylinder to gravitational radi- 
ation and the detection of gravitational waves [2] , [3] and [J . Neutron stars have 
attracted attention since it has been argued [5] that the crusts of neutron stars are 
\Q ■ in elastic states and since it has been established the existence of a solid crust and 

speculated the possibility of solid cores in neutron stars, [6], [7], [8]. 
There were many attempts to formulate a relativistic version of elasticity theory. 
Thereby laws of non relativistic continuum mechanics had to be reformulated in a 
relativistic way. The study of elastic media in special relativity was firstly carried 
out by Noether [9J in 1910 and by Born [10], Herglotz [IT] and Nordstrom [JJ] in 
1911. The discussion of elasticity theory in general relativity started with Synge 
[13] . De Witt PU, Rayner [15], Bennoun [H], [17], Hernandez [18] and Maugin 
In 1973 Carter and Quintana [20] developed a relativistic formulation of 
the concept of a perfectly elastic solid and constructed a quasi-Hookean perfect 
elasticity theory suitable for applications to high-pressure neutron star matter. Re- 
cently, Karlovini and Samuelsson [T] gave an important contribution to this topic, 
extending the results of Carter and Quintana (see also [21], [22] )• Other relevant 
formulations of elasticity in the framework of general relativity were given by Ki- 
jowski and Magli ([23], [21]) who presented a gauge- type theory of relativistic elastic 
media and a corresponding generalization [25) . The same authors also studied in- 
terior solutions of the Einstein field equations in elastic media ([20], [2"7]). 
The recent increasing consideration of relativistic elasticity in the literature shows 



1 Relativistic elasticity has been treated in the mid-20th century until the early seventies by 
many other authors. For further references, see, for example, [19] . and for later references see also 
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the win of recognition and importance of this topic, motivating for a detailed study 
of quantities used in this context, the elasticity difference tensor defined in pQ be- 
ing one of them. This tensor occurs in the relativistic Hadamard elasticity tensor 
and in the Euler equations for elastic matter. However, one can recognize the geo- 
metric role of the elasticity difference tensor, since, in principle, it can be used to 
understand the influence of the material metric (inheriting elastic properties) on 
the curvature of the space-time. 

Here, in section 2, general results about relativistic elasticity are presented. In sec- 
tion 3, the elasticity difference tensor is analysed and principal directions associated 
with this tensor are compared with the eigendirections of the pulled-back material 
metric. A specific orthonormal tetrad is introduced to write a general form of the 
elasticity difference tensor, which brings in Ricci rotation coefficients used in the 
1 + 3 formalism [25] and the linear particle densities. 

Finally, in section 4, we apply the results obtained to a static spherically symmetric 
space-time and an axially symmetric non-rotating space-time. The software Maple 
GRTensor was used to perform some calculations. 

2. General results 

Let (M, g) be a space-time manifold, i.e. a 4-dimensional, paracompact, Hausdorff, 
smooth manifold endowed with a Lorentz metric g of signature (—,+,+,+), U being 
a local chart around a point p S M. Suppose that U is filled with a continuum 
material. The material space X is an abstract 3-dimensional manifold, each point 
in X representing an idealized particle of the material. Moreover, the space-time 
configuration of the material is described by a mapping 

4- : U C M — > X, 

which associates to each point p of the space-time the particle p of the material 
which coincides with p at a certain time. Therefore \E' _1 (p) represents the flowline 
of the particle p. The operators push-forward , 3>* and pull-back *Sf* will be used 
to take contravariant tensors from M to X and covariant tensors from X to M, 
respectively, in the usual way. 

If {A = 1, 2, 3) is a coordinate system in X and {to 11 } (a = 0, 1, 2, 3) a coor- 
dinate system in U C M, then the configuration of the material can be described 
by the fields £ A = £ A (uJ a ). The mapping : T p M — > T^i p \X gives rise to a 

(3 x 4) matrix (the relativistic deformation gradient) whose entries are £ A = J^-. 
Assuming that the world-lines of the particles ^~ 1 (p) are timelike, the relativis- 
tic deformation gradient is required to have maximal rank and the vector fields 
u a 6 T P M, satisfying u a £f — 0, are required to be timelike and future oriented. 
The vector field u a is the velocity field of the matter and its components obey 
u a u a = -1, u a ^a = and u° > 0, [IS]. 

One needs to consider, in the material space X, a Riemannian metric t]ab ? describ- 
ing the "rest frame" space distances between particles calculated in the "locally 
relaxed state" or in the "unsheared state" of the material and often taken as the 
material metric. These approaches are presented in [23j and in [T], respectively. 



Capital Latin indices A,B,... range from 1 to 3 and denote material indices. Small Latin 
indices a,b,... take the values 0,1,2,3 and denote space-time indices. 
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Let tABC be the volume form of tjab, with ^f*€ABC — e abc — e abc dU d . The particle 
density form is uabc = t^ABC-, with n the particle density yielding the number 
of particles in a volume of X when integrated over that volume. One can define, 
see pp, a new tensor ItAB, which has uabc a s its volume form and is conformal to 
rjAB'- kAB — nsr/AB- This tensor will be taken as the material metric in X. 
The pull-back of the material metric 

k ab = **k AB = ggk AB (1) 

and the (usual) projection tensor 

Kb = 9ab + u a u b (2) 

are Riemannian metric tensors on the subspace of T p M orthogonal to u a . These 
tensors are symmetric and satisfy k ab u a = = h a bU a . 

The state of strain of the material can be measured by the relativistic strain tensor, 
according to e.g. [53], [2"T] : 

1 l _ 2 

Sab = ^{hab ~ Tjab) = ^{Kb ~ n 3 k ab ). (3) 

This tensor is also named as constant volume shear tensor (see [20], pQ). The 
material is said to be "locally relaxed" at a particular point of space-time if the 
material metric and the projection tensor agree at that point, i.e. if the strain 
tensor vanishes. 

When considering elastic matter sources in general relativity, one is confined to a 
stress-energy tensor taking the form T ab — —pg ab + 2t^tf = pu a u + p a b, where 

p a b = 2^^- — ph a {,, the energy density being written, for convenience, as p = ne, e 
being the energy per particle. 

Choosing an orthonormal tetrad {u, x,y, z} in M, with u in the direction of the 
velocity field of the matter and x, y, z spacelike vectors, satisfying the orthogonality 
conditions —u a u a — x a x a = y a y a — z a z a = 1, all other inner products being zero, 
the space-time metric can be written as 

9ab = -u a u b + h ab = -u a u b + x a x b + y a y b + z a z b . (4) 

Here we will choose the spacelike vectors of the tetrad along the eigendirections of 
k b = g ac k c b, so that 

k ab = n\ x a x b + njy a y b + n\ z a z b , (5) 

where n\ , and n\ are the (positive) eigenvalues of k b a . The linear particle densities 
Tlx, n% and 713 satisfy n = 711712^3 ■ It should be noticed that those eigenvectors 
are automatically orthogonal whenever the eigenvalues referred above are distinct. 
However, if the eigenvalues are not all distinct, the eigendirections associated to 
the same eigenvalue can (and will) be chosen orthogonal. 

It is convenient to consider the spatially projected connection D a acting on an 
arbitrary tensor field t b c - as follows: 

D a tt=hihl...hl..V d t e fZ . (6) 

Here V is the connection associated with g and one has D a h bc — 0. Another 
operator D, such that its action on the same tensor is 

D a tt. = Khl...hl...V d ty;- (7) 

is also considered. One has 

D b X a = D b X a + S\ C X C , (8) 
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for any space-time vector field X. The tensor field S a bc is the elasticity difference 
tensor introduced by Karlovini and Samuelsson in pQ. This third order tensor can 
be written as 

S°'f )c = —k am (Dbk mc + D c k m i) — D m kbc): (9) 

where k~ am is such that k~ am k m b = h b . This tensor is used by the same authors 
to write the Hadamard elasticty tensor, used to describe elasticity properties in 
space-time, and the Euler equations V&T ab = for elastic matter. 

The covariant derivative of the timelike unit vector field u can be decomposed as 
follows 

U a -b = ~U a Ub + D b U a = -U a llb + 7^<dh ab + (Jab + ^ab, (10) 

where ii a is the acceleration, o a p, the symmetric tracefree rate of shear tensor field, 
LiJ a f), the antisymmetric vorticity tensor field and 0, the expansion scalar field for 
the congruence associated with u. 



3. Properties of the Elasticity Difference Tensor 

Here we will investigate the algebraic properties of the elasticity difference tensor. 
This tensor, important when studying elasticity within the framework of general 
relativity, is related to the connection of the space-time, as shown in the previ- 
ous section. The following two properties of the elasticity difference tensor are 
straightforward: 

(i) it is symmetric in the two covariant indices, i. e. 

S a bc = S a c b', (11) 

(ii) it is a completely flowline orthogonal tensor field, i.e. 

S\ c u a = - S\ c u b = S a bc u c . (12) 

The elasticity difference tensor can be approached using the space-time connection, 
as will be shown here. 

It is a well known result that the difference between two connections V and V, 
associated with two different metrics g and g, respectively, defined on U, is the 
following (1,2) tensor: 

c™m; — r™m/ — r™ m ;, (13) 

T n ml and T n ml being the Christoffel symbols associated with those two metrics. In 
a local chart, this tensor can be written as ( 29 , [30 ) 

C"mi = \d np {9p m ;i + 9 P l;m ~ 5m/;p), (14) 

where g np is such that g np g pr = £™ and a semi-colon ; represents the covariant 
derivative with respect to g. The difference tensor C n ml can be used to write the 
difference of the Riemann and the Ricci tensors associated with the two metrics in 
the following form (see e.g. [5T]): 

n bed n bed — ~° bd;c "+" ° bc-d ~ u ic° bd ~r ° Zd° be \ 10 ) 

and 

Rbd — Rbd = —C a bd . a + C a ba . d — C a la C l bd + C a ld C l ba . (16) 
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The projection of the difference tensor orthogonally to u is defined by the expression 

KKh l c C n ml . (17) 

When the connections used to define the difference tensor are associated with the 
metrics g a b — —u a Ub + h a b and g a b = ~u a Ub + k a b, then the corresponding difference 
tensor, projected according to (fTT)) . yields © i.e. the elasticity difference tensor 
defined in the previous section. 

Under this approach, the elasticity difference tensor is the projection, orthogonal 
to u, of the difference between two connections, one associated with the space-time 
metric and the other with the metric g a b = —u a Ub + k a b, where k a b is the pull-back 
of the material metric Uab ■ 

Calculating the spatially projected versions of equation (fT5|) . using © and (fTT|) , 
yields the following expression for the difference of the Riemann tensors: 

hi h n g h* h q h [ K h b n h% h d q {R\ cd - R\ cd )] 
= - D eS f gh + D h Sf ge - S f ke S k gh + S f kh S k ge . (18) 



The spatially projection of (|16jl . the difference of the Ricci tensors, can be obtained 
analogously by equating the indices a = c in the last expression. 

Therefore, these expressions, which contain the elasticity difference tensor, give 
the difference between the Riemann and Ricci tensors associated with the metrics 
referred to above. 

Now we will obtain the tetrad components of the elasticity difference tensor. To 
do so, it is more convenient to use the following notation for the orthonormal 
tetrad: e° = (eg, e", eg, e|) = (u a , x a , y a , z a ). Tetrad indices will be represented 
by greek letters from the second half or the first half of the alphabet according to 
their variation as follows: /i,^, p... = — 3 and a, (3, 7... = 1 — 3. The Einstein 
summation convention and the notation for the symmetric part of tensors will only 
be applied to coordinate indices, unless otherwise stated. The operation of raising 
and lowering tetrad indices will be performed with rj^ = rf v = diag(—l, 1, 1, 1) 
and one has g a b = e lia e v wf v ■ 

Writing the Ricci rotation coefficients as y^vp — e^a^e^e^, the tetrad components 
of the elasticity difference tensor can be obtained using the standard relationship 

S a ffl = S a bc e«eX (19) 

the result being 

' - ( nr, ir: )-".:-'- i ir\ ■■■■ ir. : )-"..+(/>:- /r, )-,."-!■ L)„ ( 11;, )< ",<)'.' 

(20) 



+ D p {nl)e^-D l {nf S e la 8p 1 ]. 
An alternative form for the last expression is: 

s % = ^K 1 - e 7«)7 Q 7/3 + (1 - e/3ahVy + ^~t a ~ e ^)l^ a + m Pa 6° + m ia 5p ^ 
( n 2 \ 

where e 7a = J —j J and m a p — D a (h\n 2 p)e aa . 
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The Ricci rotation coefficients, when related to the quantities used in the decom- 
position (fTU|) . can be split into the set [52] : 

7o«o = u a (22) 
7oq/3 = g©<W + u a j3 - e Q/ 3 7 w 7 (23) 

7a/30 = -£a/3 7 ^ 7 (24) 
7a/37 = -AaSpy + Ap5 al - ^ (c^iV^ ~ e 1 spN 5 a + £ Q/3( 5 TV^). (25) 

The quantities A and N appear in the decomposition of the spatial commutation 
functions T a p = 7 Q 7(g — 7 Q g 7 , given in [33], where TV is a symmetric object. 

The elasticity difference tensor can be expressed using three second order symmetric 
tensors, here designated as M bc , a = 1, 2, 3, as follows: 



3 

S\ c = M hc x a + M bc y a + M bc z a = V M bc e a a . (26) 
1 2 3 t^i a 



Here we will study some properties of the three tensors M bc in order to understand 

a 

until which extent the principal directions of the pulled back material metric remain 
privileged directions of the elasticity difference tensor, i.e. of the tensors M bc , by 

a 

studying the eigenvalue-eigenvector problem for these second order tensors. 

First, we will obtain a general expression for M bc , a = 1,2,3, which depends 

explicitely on the orthonormal tetrad vectors, the Ricci rotation coefficients and 
the linear particle densities n a . This comes from the contraction of S a bc in ([9]) with 
each one of the spatial tetrad vectors, followed by the use of the relationships (O, 
© and appropriate simplifications. The final result is 

Mf) C ^ ( e am;(6^c) ~t~ ^(/)^Qc);m) "f" e a(6;c) &a &a(c&ab);m 
a 

+ 7o 

H r{-e™(e / 36e / 3 c n / 3n ) 3 ;rrl + e jb e 7C n 7 n 7 , m ) 

+ ™ 7 [(7o 7 a ~ 7a 7 o)u(be 7c) + e™( 

+ n p[(j0f3a - lal3a)U( b ep c ) + e™ (e^. (b e /3c) ~ ^/3{b^0cj;m)]}i 

where 7 ^ /3 ^ a, for one pair (/3, 7), a comma being used for partial derivatives. It 
should be noticed that this expression also contains the non-spatial Ricci rotation 
coefficients given in ((2"2"j) . (l2"3l and 



(27) 



Naturally, the expressions obtained for M bc still satisfy the conditions M bc u b = 0. 

a a 

The eigenvalue-eigenvector problem for M bc is quite difficult to solve in general. 

a 

However, one can investigate the conditions for the tetrad vectors to be eigenvectors 
of those tensors, the results being summarized in the two following theorems. 
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Intrinsic derivatives of arbitrary scalar fields <E>, as derivatives along tetrad vectors, 
will be represented A eQ and denned as: 

A e<s $ = $ m e™ 

where a comma is is used for partial derivatives. 

Theorem 1. The tetrad vector e a is an eigenvector for M iffn a remains invariant 

a 

along the two spatial tetrad vectors ep, such that (3 ^ a i.e. A e(3 (lnn Q ) = 
whenever f3 ^ a. 

The corresponding eigenvalue is A = A ea (lnn a ). 

Proof: In order to solve this eigenvector-eigenvalue equation the following algebraic 
conditions are used 

M£e b a e ac = A, (28) 

a 

M c h e\ e Pc = (29) 



and 



M c h e h a e lc = 0, (30) 



where 7 / |? / a. Using the orthogonality conditions satisfied by the tetrad 
vectors and the properties of the rotation coefficients, namely the fact that they 
are anti-symmetric on the first pair of indices, (|2"9"|) and (|30|) yield A efj (In n a ) = = 
A e (lnn a ) so that A e0 n a = = A e ^n a . On the other hand from (|2"5|) one obtains 
A = A ea (hm Q ). 

It should be noticed that A = whenever n a remains constant along e a . However 
this condition is equivalent to n a = c, with c a constant. In this case, k ab = 
c 2 e aa e ab + npzpaepb- □ 

Theorem 2. en is an eigenvector of M iff the following conditions are satisfied: 

a 

(i) A e/3 (lnn Q ) = 0, i.e. n a remains invariant along the direction of ep; 

(ii) Ja-fplnl ~ ™ 2 ] + 7 Q/ 3 7 [n 2 - n 2 ,] + jp ia [n^ - n|] = 0, where 7 ^ [3 ^ a for 
one pair [j3, 7). 

n 2 

The corresponding eigenvalue is A = — ^-A ea np + r y a /3p(—^- + !)• 

Proof: Contracting M£ep = \ep with e ac one obtains A e/3 (lnn a ) = 0. This 

a 

condition is satisfied whenever A ef3 n a = 0. The second condition results from 
M^e^e 7C = 0. And contracting M^ejg = Ae^ with ep c yields the eigenvalue A. The 

a a 

used simplifications are based on the orthogonality conditions of the tetrad vectors 
and on the properties of the rotation coefficients. 

Notice that the two conditions are satisfied if n a = tip = n 7 = c, where c is a 
constant. The consequence of this is that A = 0. In this 

C 2 VaVb + C 2 Z a Z b . □ 

The previous theorems show that strong conditions have to be imposed on n a , for 
a = 1,2,3, and the metric in order that the spatial tetrad vectors are principal 
directions of M, for a = 1,2,3. 



However, the conditions to have e a as eigenvector of M seem less restrictive then 

a 

the conditions for ep, for all values of /? ^ a, to be eigenvector of the same tensor 
M, since these involve not only intrinsic derivatives of the scalar fields but also 

a 

rotation coefficients of the metric. Furthermore, for e a to be an eigenvector of 
M only conditions on n a have to be satisfied, namely that n a remains constant 

a 

along the directions of ep for all values of (3 ^ a (in which case the eigenvalue 
corresponding to e a depends only on n a ). On the other hand, the conditions 
imposed for ep for all (3 ^ a to be eigenvectors of M also involve np for all j3 ^ a. 

a 

Next we will use the previous theorems to establish the conditions for e a , with 
a = 1, 2, 3 to be an eigenvector of the three tensors M, M, M simultaneously, the 

results being: 

(i) A e „(lnn Q ) = 0, 

(ii) A ea (lnnp) = 0, 

(hi) 7a/3 7 [n% - n|] + j aj p [n* - n%] + ip ia [n 2 p - n*] = 0, 
for all values of f3 and 7 such that ft ^ 7 7^ a. 

These conditions must be satisfied for all values of j3 ^ a. It is not easy to find 

the general solution to these equations, however one can say that, in general, the 

principal directions of the pulled back material metric k are not, in general, the 

principal directions of the three tensors M, M and M. It should be noticed that 

1 2 3 

the (mathematical) solution corresponding to n\ = n-z = 713 = const, is not an 
interesting result from the physical point of view. 

As a special case, we now consider that all eigenvalues of k a b are equal, i.e. 

m = n-2 = ?!3 = . (31) 

Therefore, k a b — n^h a \ )1 so that these tensors are conformally related. In physical 
terms, this corresponds to the unsheared state described in [Tj: the energy per 
particle, e, has a minimum under variations of g AB such that n is held fixed. The 
above theorems in this section simplify significantly in this case, as can easily be 
proved using (|3ip in those theorems. For completeness, we give the expressions for 
the elasticity difference tensor and the tensors M in this special case: 

a 

S\ c = ii (8 a c D b n + 5%D c n - h ad h bc D d n) , 

hhc = ^ {e ac n.b + e ab n c + (e ac u b + e ab u c )A e n - h bc A Ca n) . 
a An 

4. Examples 

Here, examples concerning the static spherically symmetric case and an axially 
symmetric, non-rotating metric are presented, where we apply the analysis devel- 
oped in the last section. The main problem when dealing with examples lies in the 
difficulties of finding an orthonormal tetrad for the space-time metric such that the 
corresponding spacelike vectors are precisely the principal directions of the pulled 
back material metric. However, in the examples presented, this difficulty was over- 
come. 
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4.1. The static spherically symmetric case. In this section we analyse the 
elasticity difference tensor and corresponding cigendirections for the static spher- 
ically symmetric metric, due to its significance on modelling neutron stars. The 
metric regarded here can be thought of as the interior metric of a non rotating star 
composed by an elastic material. 

For a static spherically symmetric spacetimc the line-clement can be written as 

where the coordinates ui a — {t, r, 9, 0} are, respectively, the time coordinate, the 
radial coordinate, the axial coordinate and the azimuthal coordinate. Choosing the 
basis one-forms u a = (-e"( r ), 0,0, 0), x a = (0, e x{ - r \ 0, 0), y a = (0,0, r, 0) and z a = 
(0, 0, 0, rsinfl) for the orthonormal tetrad, the metric is given by 

9ab — Ha'U'b ~r 

x a x b + y a yb + z a z b and h ab = x a x b + y a y b + z a z b defines the corresponding projection 
tensor. Using this tetrad, the pulled-back material metric becomes 

k ab = nfx a Xb + n\y a y b + n\z a z bl (33) 

where we have chosen n 3 — ni since for this material distribution k has only two 
different eigenvalues. 

Let l; A = {f,9,4>} be the coordinate system in the material space X. Since the 
space-time is static and spherically symmetric, f can only depend on r and one 
can take 9 — 9 and <fi = 4> so that the configuration of the material is entirely 
described by the material radius f(r). Moreover, the only non-zero components of 
the deformation gradient are -jjr = = 1 and = 1. 

In X the material metric is k^B = xaxb + TJaVb + zazb, with xa — e x drA, 
i/A = rd9 a and za — rsin9d<j)A, and where A = A(f). The pull-back of the material 
metric is then 

k\ = g ac k cb = g ac (g^k CB ) = (f ) e 2X - 2A ^ + p%5 2 b + p a 3 Sl (34) 

Comparing ([33|) and (f34|) it is simple to obtain the following values for the linear 
particle densities (all positive), which arc found to depend on r only: 

ni=ni(r) = Je X ~ A (35) 
dr 

v 

n 2 = n 2 (r) = n 3 (r) = - (36) 
r 

The non-zero components of the strain tensor ([3]), when written as functions of the 
quantities n a , are 

s rr = — e 2A (l — n^^n\) 
1 2 

s 4>4> — 7^r 2 sin 2 9(l - n~^n\) 

Using the expressions obtained for the n a one can find that the condition for this 
tensor to vanish identically is that r — ce* f Xdr . 
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Calculating the quantities given in (JTUJ) one obtains 

6 = 

0,^,0,0 
dr 

Oab'- 0-i2 — -e — — a 21 

2 dr 

2v dv 1 iv dv 
Wl2=e ^ + 2 6 Tr 

U>21 = —k>12> 

where the remaining components of a a t and LJ a b vanish. 

The non-zero components of the elasticity difference tensor S£ c are: 
f dni 
rr ~ ni dr 
g _ 1 dn 2 
n 2 dr 
f dn 2 
^ ~ n 2 dr 

C r -2A -2A n 2 -2A 2 n 2 dn 2 

S ee = re -re — -e r — — 
n\ n\ dr 

or -2A ■ 2 -2A • 2a n 2 -2A 2 ■ 2n n 2 ^ n 2 

6 ij, = e rsm — e rsm a—^ — e r sin 



Since 5' a bc = S a ci) , there are only seven non-zero components for this tensor on the 

coordinate system chosen above. 

Again, using (|35p and l|36pone obtains that: 

(i) the components S s Sr and S ^ are zero whenever the function f is of the form 
f = ar, where c± is a constant; 

(ii) S r rr is zero whenever f = c 2 + C3 j e A oV; 

(hi) the components S r es and S 17 *^ are zero whenever r = ~- . 

The second order symmetric tensors M, for a = 1, 2, 3 have the following non-zero 

a 

components: 

e A dni 

M rr = — 

1 m dr 

\ \ rck j o no, dn 2 

Mgg = e~ x r - e~ A r-f - e" V 



1 



n\ n\ dr 



M —A -2/) -A • 2/)"2 -A 2 ■ 2/)™2 <2n 2 

nf dr 



1 n i 

r <^ n 2 

M r9 = M er = — 

2 2 n 2 ar 

rsinO dn 2 
M r( h = M«A r = — 

3 3 n 2 dr 

The eigenvalues and eigenvectors of these tensors are presented in tables 1, 2 and 
3, being then compared with the eigendirections of the material metric. 
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Tabic 1 - Eigenvectors and eigenvalues for M 



Eigenvectors 




Ei 


genvalucs 




X 

y 


M2 = 


Mi 

e~ x _ 
r 


_ e~ A dm 
rti dr 

— A „2 

— % - e 

r nf 

1 - e 


-A «2 dn 2 
n\ dr 


z 


M3 = 


e~ x _ 
r 


-A n 2 dn 2 
dr 



Notice that, in the present example, M maintains the eigenvectors of k, namely x, 

y and z, the two last ones being associated with the same eigenvalue. Therefore 

the canonical form for M is Mb c — [i\XbX c + [ii{ybyc + ZbZ c ), where fi\ and p,i are 
l i 

the eigenvalues corresponding to x and y (= z), respectively. 



Table 2 - Eigenvectors and eigenvalues for M 



Eigenvectors 


Eigenvalues 


x + y 

x-y 

z 


= eZ±dn2 
~ 712 dr 

tic — _.£li*i2. 

n 2 dr 

Me = 



In this case, only the eigenvector z of k remains as eigenvector, however the corre- 
sponding eigenvalue being zero. The other two eigenvectors are x + y and x — y so 
that the canonical form for M can be expressed as M\, c = 2[i i {x}>y c + ybX c ), where 

2 2 



Table 3 - Eigenvectors and eigenvalues of M 



Eigenvectors 


Eigenvalues 


x + z 
x — z 

y 


LL-7 — £li*i2. 

M8 - n2 dr 
fig = 



Comparing M and M, it is easy to see that the role of z and y is interchanged. The 
eigenvalues of M are equal to the eigenvalues of M and the canonical form of this 
tensor field can be written as Mb c = 2^(xi } z c + ZbX c ), where ^7 = e _A (^^p — 7)- 

It should be noticed that the case ni constant is not interesting to analyze, since 
this corresponds to the vanishing of the tensors M and M. 

x, y and z would only remain eigenvectors for M and M if f would be of the form: 

r = cr, in which case M and M were reduced to a zero tensor. 
2 3 
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The tetrad components of the elasticity difference tensor can directly be obtained 
from (ED: 



_ A 1 dn\ 
n\ dr 
_> _ A 1 dn 2 



J 11 

S 2i = e 



ri2 dr 

S\ 1= e ' 



3 _ -A 1 d «2 



ri2 



J 22 



D 33 



A 1 -A 1 "1 


_ X U2 


dn 2 


2 

r r 


2 
n l 


dr 


A 1 r -\ ln l 


A n2 


dn 2 


2 

r r nf 


2 
n l 


dr 



The expressions for the Ricci rotation coefficients are 

7122 = 

r 

7133 = 

r 

cos# 

7233 = — t~z- 
rsma 



4.2. The axially symmetric non-rotating case. First, consider an elastic, ax- 
ially symmetric, uniformly rotating body in interaction with its gravitational field. 
The exterior of the body may be described by the following metric, [27], 

ds 2 = -e 2u dt 2 + e 2 >*dr 2 + e 2 "dz 2 + e 2 ^{d<P - uodt) 2 , (37) 

where v, ip, u>, fi are scalar fields depending on r and z. 

Assume that the material metric is flat. Introducing in X cylindrical coordinates 
£ A = {R, £, $}, then the material metric takes the form: 

ds 2 = dR 2 + d( 2 + R 2 d^> 2 , (38) 

where the parameters R, C depend on r and z, $ being <!>(£, r, z, <$) = <p — 0(r, z)t. 

Now, consider the limiting case of an axially symmetric non-rotating elastic system 
for which the space-time metric is given by 

This metric is obtained from (|37[) . when u> = and the angular velocity 17 = 0. 

Imposing R = R(r), £ = z and g a b — g a k(r), one obtains a further reduction to 
cylindrical symmetry. This reduction is considered in |27j . 

So, the space-time metric we will work with is given by (j39[) . where v, /i, ip depend 
on r only, and it can be written as g a b = —UaMb + % a %b + y a yb + z a Zb, where 
u a = (-e^), 0,0,0), x a = (0,e",0,0), y a = (0, 0, e"M 0) and z a = (0,0,0, e*«). 
The space-time coordinates are u a — {t, r, z, <j)}. 

In X the material metric Uab is given by Rab = xaxb + VaVb + zazb, where 
xa = dRA, Va = dzA and za — Rd<pA- The relativistic deformation gradient has 
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the following non-zero components = ^p, j~t = 1 and ^£3 = 1. Calculating 
the pull-back of the material metric one obtains 



k\ = g ac k cb = «T (£^k CB ) = e- 2 »5?5l + (^j e^S 2 + R 2 e -^6 a 3 6l (40) 
The corresponding line-element can be expressed as 

ds 2 = dr 2 +( ^f) dz 2 + R 2 d<f> 2 . (41) 



On the other hand, the material metric in the space-time M is given by 

k a b = n\x a x b + n%y a y b + n\z a z h . (42) 

Comparing (|4T)|) with (|4"2"]) one concludes that the linear particle densities (all posi- 
tive) are expressed as 



rix =ni(r) = e-* (43) 

n 2 =n 2 (r) = e - ' 1 -? (44) 
dr 

n 3 = n 3 (r) = Re~^ . (45) 



The strain tensor ([3]) is composed of the following components 

1 

— ( 

2 

1 



s rr = -e 2fJ -(l - n s/if) 



s zz = -e 2 »(l- n -*nl) 
S0 = -e 2 ^{l -n~in 2 ) 



The strain tensor vanishes if the condition R(r) = r = e^-f 1 is satisfied. 
Calculating the quantities given in (jTUJ) one obtains 

9 = 

Ua= (°' e2l/ ^' '° 

&ab- fl2 = -e — - = 021 

2 ctr 

2v dv 1 4 „diy 
Wl2 =e - + -e - 

where the remaining components of a ab and uj ab vanish. 
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The non-zero components of the elasticity difference tensor are 

_ 1 dni 
rr ~ m dr 
1 dn 2 
zr n 2 dr 
c<t> 1 dn 3 
^ ~ n 3 dr 
gr _ dfi n\ dfx n 2 dn 2 
zz dr n\ dr n\ dr 



or 

O A 



e 



-20-2^ fdjp _ nicty _ n 3 dn 3 



dr n\ dr n\ dr 



It can be observed that only seven components of the elasticity difference tensor 
are non-zero. 

Using the expressions (|43j). (|44p and (|45jl one can conclude that: 

(i) S r rr is zero whenever /i(r) = c, where c is a constant; 

(ii) S z zr is zero whenever R(r) = c\ + c 2 J e^ r 'dr; 

(iii) S"^ r is zero whenever R(r) — c 3 e^^\ 

(iv) S r zz is zero whenever R(r) = ± J y/2fi(r) + c<idr + C5; 

( v ) s ' r 00 is zer0 whenever R(r) = ±^2 J ^r^dr + c 6 . 

The second-order tensors M, M and M have the following non-zero components: 

12 3 



M rr 
1 


1 dm 

ni dr 




M zz 
1 


_ efl ( dfi n\ 

\ dr At 


n 2 

9 
nf 


1 


\ dr n\ 


di/> 
dr 


M rz 
2 


2 n 2 dr 




M rt j, 
3 


= M, r = e^-^ 
3 "3 dr 





n\ dr 



The next three tables contain the eigenvalues and eigenvectors for these tensors, 
which are then compared with the eigenvectors of the pulled-back material metric. 



Table 1 - Eigenvectors and eigenvalues for M 



Eigenvectors 


Eigenvalues 


X 




- 1 n± dr 




V 


A 2 


V dr dr 


n 2 dn-2 \ 
dr J 


z 


A 3 


— e -p, (£>£_ ntdjj, _ 
\ dr dr 


n 3 dn 3 \ 
dr J 



One can observe that the eigendirections x, y and z of k are also eigenvectors for 
the tensor M and the eigenvectors are associated with different eigenvalues. The 

canonical form for M can be written as Mt, c = \\XhX c + X 2 ybVc + X 3 ZbZ c . 
1 1 
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Table 2 - Eigenvectors and eigenvalues for M 



Eigenvectors 


Eigenvalues 


x + y 

x-y 

z 


A 4 = e"" 

A 5 = -e-" 
A 6 = 


f 1 dn 2 \ 
1 n 2 rfr y 

/ _J_ ring 
\ 7i2 dr 

-- 


) 



M inherits only the eigenvector z of k, which corresponds to a zero eigenvalue. The 
other two eigenvectors of M are linear combinations of x and y: x + y and x — y, 
whose corresponding eigenvalues are symmetric in sign. The canonical form for M 

can be written as M bc = 2\ i (x b y c + ybX c ), where A 4 = I -fe- - -£; 



Table 3 - Eigenvectors and eigenvalues for M 



Eigenvectors 


Eigenvalues 


x + z 
x — z 

y 


A 7 = e"" 

A 8 = -e-" 
A 9 = 


f 1 drt 3 \ 
^n 3 cir y 

/ 1 rfn 3 
y n 3 dr 

= 


) 



M inherits the eigenvalue y of k, which is associated with the eigenvalue zero. 
The other two eigenvectors of M are linear combinations of x and z: x + z and 

x — z. These two eigenvectors are associated with sign symmetric eigenvalues. 

The canonical form for M can be written as Mf, c — 2Xy(xi,z c + Zf,x c ), where A7 = 
3 3 



( L^R _ dtp \ 

\R dr dr J 



x and y would only be eigenvectors for M if R(r) would be of the form R(r) — 

c\ + J e^drc2, but in this case M would vanish, x and z would only be eigenvectors 

for M if R(r) would be of the form R(r) = c^e? and this would reduce M to a zero 
tensor. 

One can observe that the role that y and n 2 play for the tensor M is the same that 
z and 713 play for M. That is, the results for M and M are very similar, only y 
and n 2 are substituted by z and ^13, respectively. 



The tetrad components of the elasticity difference tensor obtained from (|2Tj) and 
the expressions for the Ricci rotation coefficients are listed below: 
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= e 


-fj, 1 dn i 
n\ dr 










c2 
° 21 


= e 


n 2 dr 










n3 
° 31 


= e 


- M 1 dn 3 
TI3 dr 










° 22 


= e 


dr 


-M n 2 

n? 


dr 


2 


dn 2 
dr 


° 33 


= e 


dr 


9 
Hi 




9 

nf 


dn 3 
dr 
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The expressions for the Ricci coefficients are 



a, dr 

dr 

7133 - — • 
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